NMR structure of the thrombin-binding DNA aptamer stabilized by Sr2+.

The structure of thrombin-binding DNA aptamer complexed with a single Sr2+ ion (Sr2+:TBA complex) has been determined using NMR spectroscopy and restrained molecular dynamics simulations. The quadruplex structure for the Sr2+:TBA complex is similar in topology, but distinct in structure, from that previously reported for the K+:TBA complex. The inter-tetrad distance of the Sr2+:TBA complex is 3.8 angstroms, or 0.7 angstroms larger than in the K+:TBA complex. This substantial difference can be attributed to a different binding site for Sr2+ in the Sr2+:TBA complex than for K+ in the K+:TBA complex. The Sr2+:TBA complex assumes a 1:1 stoichiometry, and it is very likely that the Sr2+ ion simultaneously interacts with the eight O6 atoms of the two G-tetrads. The results indicate that quadruplex DNA structures are highly sensitive to the presence of specific metal ions. The binding of specific metal ions may modulate the biological activity of quadruplex DNA structures in vivo.

Folding of the thrombin aptamer into a G-quadruplex with Sr(2+): stability, heat, and hydration.

It has been shown that the DNA aptamer d(G(2)T(2)G(2)TGTG(2)T(2)G(2)) adopts an intramolecular G-quadruplex structure in the presence of K+. Its affinity for trombin has been associated with the inhibition of thrombin-catalyzed fibrin clot formation. In this work, we used a combination of spectroscopy, calorimetry, density, and ultrasound techniques to determine the spectral characteristics, thermodynamics, and hydration effects for the formation of G-quadruplexes with a variety of monovalent and divalent metal ions. The formation of cation-aptamer complexes is relatively fast and highly reproducible. The comparison of their CD spectra and melting profiles as a function of strand concentration shows that K+, Rb+, NH(4)+, Sr(2+), and Ba(2+) form intramolecular cation-aptamer complexes with transition temperatures above 25 degrees C. However, the cations Li+, Na+, Cs+, Mg(2+), and Ca(2+) form weaker complexes at very low temperatures. This is consistent with the observation that metal ions with ionic radii in the range 1.3-1.5 A fit well within the two G-quartets of the complex, while the other cations cannot. The comparison of thermodynamic unfolding profiles of the Sr(2+)-aptamer and K+ -aptamer complexes shows that the Sr(2+)-aptamer complex is more stable, by approximately 18 degrees C, and unfolds with a lower endothermic heat of 8.3 kcal/mol. This is in excellent agreement with the exothermic heats of -16.8 kcal/mol and -25.7 kcal/mol for the binding of Sr(2+) and K+ to the aptamer, respectively. Furthermore, volume and compressibility parameters of cation binding show hydration effects resulting mainly from two contributions: the dehydration of both cation and guanine atomic groups and water uptake upon the folding of a single-strand into a G- quadruplex structure.

Incorporation of a cationic aminopropyl chain in DNA hairpins: thermodynamics and hydration.

We report on the physicochemical effects resulting from incorporating a 5-(3-aminopropyl) side chain onto a 2'-deoxyuridine (dU) residue in a short DNA hairpin. A combination of spectroscopy, calorimetry, density and ultrasound techniques were used to investigate both the helix-coil transition of a set of hairpins with the following sequence: d(GCGACTTTTTGNCGC) [N = dU, deoxythymidine (dT) or 5-(3-aminopropyl)-2'-deoxyuridine (dU*)], and the interaction of each hairpin with Mg(2+). All three molecules undergo two-state transitions with melting temperatures (T(M)) independent of strand concentration that indicates their intramolecular hairpin formation. The unfolding of each hairpin takes place with similar T(M) values of 64-66 degrees C and similar thermodynamic profiles. The unfavorable unfolding free energies of 6.4-6.9 kcal/mol result from the typical compensation of unfavorable enthalpies, 36-39 kcal/mol, and favorable entropies of approximately 110 cal/mol. Furthermore, the stability of each hairpin increases as the salt concentration increases, the T(M)-dependence on salt yielded slopes of 2.3-2.9 degrees C, which correspond to counterion releases of 0.53 (dU and dT) and 0.44 (dU*) moles of Na(+) per mole of hairpin. Absolute volumetric and compressibility measurements reveal that all three hairpins have similar hydration levels. The electrostatic interaction of Mg(2+) with each hairpin yielded binding affinities in the order: dU > dT > dU*, and a similar release of 2-4 electrostricted water molecules. The main result is that the incorporation of the cationic 3-aminopropyl side chain in the major groove of the hairpin stem neutralizes some local negative charges yielding a hairpin molecule with lower charge density.

Shape-selective binding of geometrically-constrained bis-distamycins to a DNA duplex and a model Okazaki fragment of identical sequence.

The binding of ligands to nucleic acids is of great interest for the control of gene expression and other nucleic acid mediated processes. We have evaluated the binding of several geometrically-constrained bis-distamycins to a model Okazaki fragment [OKA], or a DNA duplex having identical base sequence [DD], using gel-shift assays, optical spectroscopy and differential scanning calorimetry. In the case of covalent attachment of two distamycins to a central benzene ring, a similar binding profile was observed for [DD] as was observed for [OKA] (para binds [K(app) > 10(6) M(-1)], meta binds only weakly). For a central pyridyl ring, however, clear distinction between the binding to [DD] and binding to [OKA] was observed. While none of the three meta isomers having a central pyridyl ring bound [OKA], two of them (MT-17 and MT-12) bound [DD] [K(app) > 10(6) M(-1)]. These results demonstrate subtle differences in lexitropsin shape and placement of electronegative atoms may result in selective binding to a nucleic acid duplex based both on base sequence and chemical composition. Selective binding to DNA duplexes may be useful for designing ligands that regulate transcription, but do not interfere in other nucleic acid mediated processes.

Shape-selective recognition of a model Okazaki fragment by geometrically-constrained bis-distamycins.

Okazaki fragments represent interesting targets for the design of anticancer drugs because of their selective occurrence during DNA replication, a process often elevated in aggressive malignancies. Structural studies have indicated a bend occurs in the helical axis at the junction region (JR) that joins the DNA duplex region (DDR) and the RNA-DNA hybrid duplex region (HDR) of model Okazaki fragments. To identify a structural motif that provides a shape complementary to the Okazaki fragment minor groove, we have investigated the binding of geometrically-constrained bis-distamycins to a model Okazaki fragment, [OKA], with a sequence derived from the genome of simian virus 40 (SV40). Both the JR and the DDR of [OKA] contain consecutive A/T base pairs that could accommodate distamycin binding. Of the six bis-distamycins selected for analysis, the two with a para configuration of the distamycins on the benzene or pyridine scaffold bound [OKA] tightly (Kd approximately 10(-6) M from gel-shift assays; Kd approximately 10(-8) M from deltaT(M)) while the four with a meta orientation did not bind. The two mono-distamycins studied also did not bind [OKA]. Molecular modeling of the complex between the para bis-distamycin MT-9 and [OKA] revealed MT-9 adopted an S- shape complementary to the minor groove of the model Okazaki fragment.

[Volume effect and change in compressibility in intercalation of ethidium bromide and daunomycin in DNA]. / Ob''emnyi éffekt i izmenenie szhimaemosti pri interkaliatsii étidiuma bromida i daunomitsina v DNK.

A combination of high-precision and density techniques was used to investigate the binding of ethidium bromide and daunomycin to DNA. It was shown that intercalation is accompanied by a decrease in apparent volume and adiabatic compressibility of DNA, which were -28 cm3mol-1 and -20.10(-4) cm3mol-1bar-1 for binding of ethidium to DNA, and -28 cm3mol-1 and -18.10(-4) cm3mol-1bar-1 for daunomycin, respectively. The values obtained have been compared with those of the hydration contributions of atomic groups to the apparent molar volume and apparent molar compressibility of the ligands as well as the contribution of intrinsic compressibility of DNA. The volume effects in the intercalation of ligands correspond to the values expected for the dehydration of aromatic heterocyclic rings. On other hand, negative compressibility changes may be explained by a decrease of the intrinsic compressibility of DNA and possible increase of its hydration as result of conformational changes upon ligand binding.

Thermodynamics of dT--dT base pair mismatching in linear DNA duplexes and three-arm DNA junctions.

We have used a combination of magnetic-suspension densimetry and calorimetry to derive complete thermodynamic profiles, including volume changes, for the formation of linear DNA duplexes and three-arm branched DNA junctions, from their component strands, with and without dT-dT mismatches. The formation of each type of complex at 20 degrees C is accompanied by a favorable free energy, with a favorable enthalpy term partially compensated by an unfavorable entropy. Formation is associated also with net uptake of water molecules. Using the formation of the fully-paired linear duplex or three-arm junction as reference states, we can establish a thermodynamic cycle in which the contribution of the single-strand species cancels. From this cycle, we determine that substitution of dA for dT has a differential free energy of deltadeltaG degrees of +2.4 kcal mol(-1) for mismatched duplex and +2.0 kcal mol(-1) (on the average) for the mismatched junction. These unfavorable differential free energies result from an unfavorable enthalpy, partially compensated by a favorable entropy, and a negative deltadeltaV. The free energies in the two cases have signs opposed to those of deltadeltaV, a situation that implicates hydration changes in creating the mismatch. When the deltadeltaV terms are normalized by the total number of base pairs involved, the immobilization of structural water molecules (and/or substitution of electrostricted for hydrophobic water molecules) is about 7 times greater for junctions than duplexes. This is consistent with more extensive hydrophobic hydration of branched DNA structures than of duplexes.

Thermodynamic investigation of the association of ethidium, propidium and bis-ethidium to DNA hairpins.

We have used a combination of calorimetric and spectroscopic techniques to investigate the association of the bis-intercalator ethidium homodimer (bis-ethidium) to short DNA hairpins with sequences: d(GCGCT5GCGC) and d(CGCGT5CGCG). The helix-coil transition of each hairpin, investigated by UV and calorimetric melting protocol, takes place in monomolecular two-state transitions with characteristic enthalpies of approximately 37 kcal mol-1 for disrupting the four dG-dC base pairs of the hairpin stems. Deconvolution of the bis-ethidium-hairpin calorimetric titration curves indicate that each hairpin contains two distinct binding sites for the ligand: a high affinity site in the stem (Kb approximately 10(7)) that accommodates one bis-ethidium molecule and a lower affinity site (Kb approximately 10(6)) located probably at the loop that accommodates two bis-ethidium molecules. The overall stoichiometries of three ligands per hairpin are in agreement with those obtained in continuous variation experiments using visible spectroscopy. The interaction of bis-ethidium for each type of sites results in enthalpy driven reactions, with average binding enthalpies, delta Hb, of -13.1 and -12.1 kcal mol-1 for the stem and loop sites, respectively. Comparison to the thermodynamic profiles of ethidium and propidium binding reveals that the bis-ethidium binding to the stem site of each hairpin has a more favorable free energy term of -1.4 kcal mol-1 and more favorable enthalpy of -4.2 kcal mol-1. These suggest that only one phenanthridine ring of bis-ethidium intercalates in the stem, while the second planar ring is exposed to solvent or weakly associated to the surface of DNA.

Interaction of minor groove ligands to an AAATT/AATTT site: correlation of thermodynamic characterization and solution structure.

A combination of circular dichroism spectroscopy, titration calorimetry, and optical melting has been used to investigate the association of the minor groove ligands netropsin and distamycin to the central A3T2 binding site of the DNA duplex d(CGCAAATTGGC).d(GCCAATTTGCG). For the complex with netropsin at 20 degrees C, a ligand/duplex stoichiometry of 1:1 was obtained with Kb approximately 4.3 x 10(7) M-1, delta Hb approximately -7.5 kcal mol-1, delta Sb approximately 9.3 cal K-1 mol-1, and delta Cp approximately 0. Previous NMR studies characterized the distamycin complex with A3T2 at saturation as a dimeric side-by-side complex. Consistent with this result, we found a ligand/duplex stoichiometry of 2:1. In the current study, the relative thermodynamic contributions of the two distamycin ligands in the formation of this side-by-side complex (2:1 Dst.A3T2) were evaluated and compared with the thermodynamic characteristics of netropsin binding. The association of the first distamycin molecule of the 2:1 Dst.A3T2 complex yielded the following thermodynamic profile: Kb approximately 3.1 x 10(7) M-1, delta Hb = -12.3 kcal mol-1, delta Sb = -8 cal K-1 mol-1, and delta Cp = -42 cal K-1 mol-1. The binding of the second distamycin molecule occurs with a lower Kb of approximately 3.3 x 10(6) M-1, a more favorable delta Hb of -18.8 kcal mol-1, a more unfavorable delta Sb of -34 cal K-1 mol-1, and a higher delta Cp of -196 cal K-1 mol-1. The latter term indicates an ordering of electrostricted and structural water molecules by the complexes. These results correlate well with the NMR titrations and are discussed in context of the solution structure of the 2:1 Dst.A3T2 complex.

Equilibrium analysis of ethidium binding to DNA containing base mismatches and branches.

In the processes of DNA replication, recombination, and repair, duplex DNA can transiently form branched structures, such as Holliday junctions, as well as base pair mismatches and bulges. These stages have altered ligand and protein binding properties from normal double helical DNA. A variety of ligands have been reported to interact more tightly at branches and bulges than to normal duplex sites. The stoichiometry, structural basis, and thermodynamics of this effect have not been determined. We have investigated the binding of the intercalator, ethidium bromide, to several DNA constructs including base mismatches, bulges, and three- and four-arm branched structures, using chemical footprinting, titration calorimetry, and fluorescence lifetime measurements. Two classes of binding sites are detected in three- and four-arm junctions in our high ionic strength conditions: one class is characterized by a small number of ligands (2-4 per DNA), with high binding affinity (K > 10(5)), and the second by a larger number of sites (10-12 per DNA) with lower affinity (K approximately 10(4)). By use of appropriate control experiments, the former appear to be associated with sites at or near the branch point or mismatch, while the latter are consistent with binding to the normal duplex DNA region(s) of the molecule. Titration calorimetry indicates an enthalpy of -10 to -13 kcal/mol for binding of ethidium to a mismatch or three- and four-arm branch point. The tight binding class is associated with a fluorescence lifetime of 12-16 ns, distinct from that of free ethidium (ca. 2 ns) and the longer lifetime observed for ethidium intercalated in duplex DNA (22-26 ns).

Differential hydration of dA.dT base pairs in parallel-stranded DNA relative to antiparallel DNA.

Parallel-stranded DNA is a novel double-stranded helical form of DNA. Its secondary structure is established by reverse Watson-Crick base pairing between the bases of the complementary strands forming a double helix with equivalent grooves. We have used a combination of magnetic suspension densitometry and isothermal titration calorimetry to obtain complete thermodynamic profiles for the formation of two DNA 25 mer duplexes. The duplexes contain exclusively dA.dT base pairs in either parallel (ps-d1.D2) or antiparallel (aps-D1.D3) orientation. At 15 degrees C, the formation of each duplex is accompanied by favorable free-energy terms resulting from the partial compensation of favorable exothermic enthalpies and unfavorable entropies and by an uptake of both counterions and water molecules. By taking into account the contribution of single-strand base-stacking interactions and using the formation of the aps-D1.D3 duplex as a reference state to establish a thermodynamic cycle in which the similar single strands cancel out, we obtained a delta delta G zero term of +10 kcal mol-1 duplex formed that results from a partial differential enthalpy-entropy compensation of +32 kcal mol-1 and a delta delta V of 257 mL mol-1. The positive sign of this enthalpy-entropy compensation together with the marginal differential counterion uptake of 0.2 mol of Na+/mol of duplex is characteristic of processes driven by differential hydration and strongly suggests that the parallel duplex is much less hydrated than its antiparallel counterpart by approximately 4 mol of water/mol of base pair.

Thermodynamics of DNA hairpins: contribution of loop size to hairpin stability and ethidium binding.

A combination of calorimetric and spectroscopic techniques was used to evaluate the thermodynamic behavior of a set of DNA hairpins with the sequence d(GCGCTnGCGC), where n = 3, 5 and 7, and the interaction of each hairpin with ethidium. All three hairpins melt in two-state monomolecular transitions, with tm's ranging from 79.1 degrees C (T3) to 57.5 degrees C (T7), and transition enthalpies of approximately 38.5 kcal mol-1. Standard thermodynamic profiles at 20 degrees C reveal that the lower stability of the T5 and T7 hairpins corresponds to a delta G degree term of +0.5 kcal mol-1 per thymine residue, due to the entropic ordering of the thymine loops and uptake of counterions. Deconvolution of the ethidium-hairpin calorimetric titration curves indicate two sets of binding sites that correspond to one ligand in the stem with binding affinity, Kb, of approximately 1.8 x 10(6) M-1, and two ligands in the loops with Kb of approximately 4.3 x 10(4) M-1. However, the binding enthalpy, delta Hb, ranges from -8.6 (T3) to -11.6 kcal mol-1 (T7) for the stem site, and -6.6 (T3) to -12.7 kcal mol-1 (T7) for the loop site. Relative to the T3 hairpin, we obtained an overall thermodynamic contribution (per dT residue) of delta delta Hb = delta(T delta Sb) = -0.7(5) kcal mol-1 for the stem sites and delta delta Hb = delta(T delta Sb) = -1.5 kcal mol-1 for the loop sites. Therefore, the induced structural perturbations of ethidium binding results in a differential compensation of favorable stacking interactions with the unfavorable ordering of the ligands.

Contribution of loops and nicks to the formation of DNA dumbbells: melting behavior and ligand binding.

We have evaluated the thermodynamic contribution of thymine loops and nicks to the overall stability of double-helical DNA by investigating (1) the melting behavior of two unligated DNA dumbbells and their corresponding core duplexes and (2) the association of netropsin to the central core of four A.T base pairs of these molecules. Temperature-dependent UV absorption and differential scanning calorimetry techniques have been used to characterize the helix-coil transitions of all four deoxyoligonucleotide duplexes. In 10 mM NaP(i) buffer at pH 7.0, all transitions were monophasic. The dumbbells melt with transition temperatures, Tm, independent of strand concentration, while each duplex melts with transition temperature dependence on strand concentration, characteristic of mono- and bimolecular processes, respectively. The Tm's for the dumbbells correspond to those of single hairpins containing only four base pairs in the stem. We obtain dTm/d log [Na+] values of 10.9-12.5 degrees C for these molecules, which correspond to similar counterion releases and suggest helical structures with similar charge densities and helical strandedness. Standard thermodynamics profiles at 5 degrees C reveal that the favorable free energy of forming these ordered structures results from the partial compensation of favorable enthalpies with unfavorable entropies. The stabilization of the dumbbells relative to the core duplexes is enthalpic, due to extra stacking of the nearest loop thymines on the G.C base pairs at both ends of the stem.(ABSTRACT TRUNCATED AT 250 WORDS)

Volume changes correlate with entropies and enthalpies in the formation of nucleic acid homoduplexes: differential hydration of A and B conformations.

We have used a combination of densimetric, calorimetric, and uv absorption techniques to obtain a complete thermodynamic characterization for the formation of nucleic acid homoduplexes of known sequence and conformation. The volume change delta V accompanying the formation of four duplexes was interpreted to reflect changes in hydration based on the electrostriction phenomenon. In 10 mM sodium phosphate buffer at pH 7, the magnitude of the measured delta V's ranged from -2.0 to +7.2 ml/mol base pair and followed the order of poly(rA).poly(dT) approximately poly(dA).poly(dT) < poly(rA).poly(dU) approximately poly(rA).poly(rU). Inclusion of 100 mM NaCl in the same buffer gave the range of -17.4 to -2.3 mL/mol base pair and the following order: poly(dA).poly(dT) < poly(rA).poly(dT) < poly(rA).poly(rU) approximately poly(rA).polyr(dU). Standard thermodynamic profiles of forming these duplexes from their corresponding complementary single strands indicated similar free energies that resulted from the compensation of favorable enthalpies with unfavorable entropies along with a similar counterion uptake at both ionic strengths. The differences in these compensating effects of entropy and enthalpy correlated very well with the volume change measurements in a manner suggesting that the homoduplexes in the B conformation are more hydrated than are those in the A conformation. Moreover, the increased thermal stability of these homoduplexes resulted from an increase in the salt concentration corresponding to larger hydration levels as reflected by the delta V results.

T-T base mismatches enhance drug binding at the branch site in a four-arm DNA junction.

Base mismatches--non Watson-Crick pairing between bases--can arise in duplex DNA as a consequence of mutational events or by recombination. In a duplex, the sequence of the two bases involved, and those flanking the site of mismatch, determines the local structure and extent of destabilization of the helix. Base mismatches can arise also in recombination of nonhomologous strands, and their occurrence in Holliday recombination intermediates can influence the outcome of general or specialized recombination events. We have previously reported that the branch site in a DNA junction can interact selectively with a variety of ligands. Here we describe the thermodynamics of junctions containing T-T mismatches flanking the branch and show that these structures bind methidium and other intercalators with higher affinity than junctions lacking mismatches.

Differential hydration of homopurine sequences relative to alternating purine/pyrimidine sequences.

The minor groove ligand distamycin A has been used to probe the relative hydration of the minor groove of eight synthetic polynucleotides of known sequence and composition. A combination of densimetric, calorimetric, and temperature-dependent spectroscopic techniques have been used to obtain complete thermodynamic profiles (delta Gzero, delta Hzero, delta Szero, and delta Vzero) for the association of distamycin A to all polymer duplexes. In 10 mM phosphate buffer, pH 7, binding of the drug to each of the polymeric duplexes resulted in characteristic negative changes in both the volume and enthalpy. Although the binding constants were found to be identical for pairs of isomer polynucleotides having identical compositions but different sequences, the values of delta Hzero, delta Szero, and delta Vzero of each such pair were remarkably different. The entropy changes were found to roughly parallel the volume changes; no such trend was seen between delta Hzero and delta Vzero. The data support the hypothesis that the volume changes observed for these systems reflect the coulombic-hydration contribution to the entropy. The heteropolymer duplexes generated much larger exothermic contributions, less favorable entropies and larger volume contractions than did the corresponding homopolymer duplexes of identical composition, and strongly suggest that polynucleotides with homopurine sequences are more hydrated than polynucleotides with alternating purine/pyrimidine sequences. In addition, it was found that duplexes containing guanine sharply reduced the affinity for the drug, also lowering the exothermicity but raising the entropy. This may be explained by the presence of an amino group in the minor groove that prevents hydrogen bonding. Substitution of the guanine with inosine reversed this trend in the thermodynamic properties. Furthermore, substitution of poly(dA) for poly(rA) in a duplex produced a similar reduction in the affinity, while raising the exothermic contribution and greatly reducing the favorable entropy effect in agreement with an apparent increase in the hydration state.

Interaction of the dye ethidium bromide with DNA containing guanine repeats.

DNA containing one or more copies of the motifs repeated in telomere sequences has unusual conformational properties. The isolated sequence from the protozoan Oxytricha, dT4G4 has the potential to form tetramers in the presence of sodium or potassium ions. We report here that these tetramers bind ethidium tightly, with an interaction that fulfills several criteria for an intercalative mechanism in the G sequence. By contrast, the 4-fold tandem repeat of this subunit, d(T4G4)4, does not interact specifically with ethidium in the presence of Na+. This difference might have a simple structural basis: the tetramer of dT4G4 forms a stack of four G-quartets in the presence of Na+ or K+, whereas the constraint imposed by the T4 "tethers" in the repeat d(T4G4)4 allows only two layers to form in the presence of Na+. In the presence of sufficient K+, the latter can partially form a four-layer G-quartet structure, which interacts with ethidium. This idea is supported by analysis of a "relaxed" sequence, dT4G4(T7G4)3, which allows formation of four G-quartets and binds ethidium in the presence of Na+ as well as K+. Ethidium (and intercalators generally) should thus be able to retard or inhibit the action of telomerase in the presence of K+.

Thermodynamics of DNA branching.

Branched DNA molecules arise transiently as intermediates in genetic recombination or on extrusion of cruciforms from covalent circular DNA duplexes that contain palindromic sequences. The free energy of these structures relative to normal DNA duplexes is of interest both physically and biologically. Oligonucleotide complexes that can form stable branched structures, DNA junctions, have made it possible to model normally unstable branched states of DNA such as Holliday recombinational intermediates. We present here an evaluation of the free energy of creating four-arm branch points in duplex DNA, using a system of two complementary junctions and four DNA duplexes formed from different combinations of the same set of eight 16-mer strands. The thermodynamics of formation of each branched structure from the matching pair of intact duplexes have been estimated in two experiments. In the first, labeled strands are allowed to partition between duplexes and junctions in a competition assay on polyacrylamide gels. In the second, the heats of forming branched or linear molecules from the component strands have been determined by titration microcalorimetry at several temperatures. Taken together these measurements allow us to determine the standard thermodynamic parameters for the process of creating a branch in an otherwise normal DNA duplex. The free energy for reacting two 16-mer duplexes to yield a four-arm junction in which the branch site is incapable of migrating is + 1.1 (+/- 0.4) kcal mol-1 (at 18 degrees C, 10 mM-Mg2+). Analysis of the distribution of duplex and tetramer products by electrophoresis confirms that the free energy difference between the four duplexes and two junctions is small at this temperature. The associated enthalpy change at 18 degrees C is +27.1 (+/- 1.3) kcal mol-1, while the entropy is +89 (+/- 30) cal K-1 mol-1. The free energy for branching is temperature dependent, with a large unfavorable enthalpy change compensated by a favorable entropy term. Since forming one four-stranded complex from two duplexes should be an entropically unfavorable process, branch formation is likely to be accompanied by significant changes in hydration and ion binding. A significant apparent delta Cp is also observed for the formation of one mole of junction, +0.97 (+/-0.05) kcal deg-1 mol-1.